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Intrinsic versus Super-rough anomalous scaling in Spontaneous Imbibition 
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We study spontaneous imbibition using a phase field model in a two dimensional system with a 
dichotomic quenched noise. By imposing a constant pressure /X;, < at the origin, we study the 
case when the interface advances at low velocities, obtaining the scaling exponents z = 3.0 ± 0.1, 
Q = 1.50 ± 0.02 and aioc = 0.95 ± 0.03 within the intrinsic anomalous scaling scenario. These 
results are in quite good agreement with experimental data recently published. Likewise, when we 
increase the interface velocity, the resulting scaling exponents are z = 4.0 ±0.1, a = 1.25 ±0.02 and 
O-ioc = 0.95 ± 0.03. Moreover, we observe that the local properties of the interface change from a 
super-rough to an intrinsic anomalous description when the contrast between the two values of the 
dichotomic noise is increased. From a linearized interface equation we can compute analytically the 
global scaling exponents which are comparable to the numerical results, introducing some properties 
of the quenched noise. 

PACS numbers: 68.35.Ct 



I. INTRODUCTION 



The growth of fluctuating interfaces in disordered me- 
dia has been a subject of much interest in last few 
years 0, 13 • The study of the interface kinetic roughen- 
ing is a fundamental problem in non-equilibrium statisti- 
cal physics and has important applications in technolog- 
ical processes and material characterization. Examples 
of experimental studies are fracture surfaces 0, Q , slow 
combustion fronts 0, fluid-air interfaces in porous me- 
dia IE S S IS The main common feature of the 
stable fluctuating interface h(x,t) in such experiments 
is that it can be described in terms of scaling laws [lll| . 
The basic assumption is to consider the Family- Vicsek 
hypothesis ^3 1 which supposes a self-afflne interface de- 
scribed by a growth exponent (3, a roughness exponent 
a and a dynamic exponent z, with the scaling relation 
a = zf3. Although this scaling hypothesis is valid in a 
great variety of cases, different experiments0,|3 and nu- 
merical models |lj,|i3 suggest a new scaling description, 
the so-called anomalous scaling, in the sense that the in- 
terface behaves in a not self-afflne manner and the 
local and global interface fluctuations are quite different 
so that a new local roughness exponent aioc has to be in- 
troduced. As pointed out by J.M. Lopez anomalous 
scaling appears due to the fact that the mean local slope 
of the interface (V/i) has a nontrivial dynamics. Further- 
more, depending on the local exponent values, anomalous 
scaling can be classified in a different groups such as the 
super-rough or intrinsic anomalous scaling |l7l |. 

The advancement of a fluid into a porous medium has 
received much attention both experimentally and numer- 
ically P, • It is important to remark that experiments 
can be carried out in two different ways: a) forced-flow 
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imbibition, where the invading fluid gets into the porous 
medium at a constant injection rate and b) spontaneous 
imbibition, where the fluid advances into the medium due 
to the capillarity. In the present paper we focus on the 
spontaneous imbibition phenomenon ( see Sec. Ill Al ), 
and we base our study on an experimental work appeared 
recently where an oil-air interface advances in a hori- 
zontal Hele-Shaw cell with a random gap spacing between 
both glass plates. Traction due to capillarity is so strong 
that an external negative pressure must be applied to 
have a very slow motion of the roughening front. With 
these experimental conditions, the scaling exponents ob- 
tained are z = 3 , a ~ 2 and aioc — 1 within the intrinsic 
anomalous scaling framework. 

The first attempt to describe spontaneous imbibition 
numerically was proposed by Dube et al pol l2ll | employ- 
ing a phase field model Hi IH HI 13 . The model is 
based on the introduction of a partial differential equa- 
tion (model B in |23|) with a quenched noise and it in- 
cludes the local conservation law of the liquid bulk mass. 
The resulting scaling exponents in the studied regime 
are z = A, a = 1.25 and aioc — 1 within the super-rough 
anomalous scaling. 

The objective of the present paper is to make a deeper 
study of the spontaneous imbibition using the phase field 
model. Our aim is to reproduce the experimental result 
13 and explain the origin of the intrinsic anomalous scal- 
ing. It turns out that the velocity of the interface is 
a relevant parameter. The exponent z = 3 indicates a 
global dynamics driven by capillary forces which occurs 
at low velocities. Intrinsic anomalous scaling has never 
numerically been observed in spontaneous imbibition and 
a dynamical growth with z = 3 has only been observed 
previously in the forced imbibition case P, ^E 113 • 
increasing the interface velocity we recover the z = 4 pre- 
dicted by Dube et al 113 • In the same line, by increas- 
ing the capillary contrast of the quenched noise we can 
move from a super-rough to an intrinsic anomalous seal- 
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ing. Furthermore, the scaling properties of the quenched 
noise will be discussed to determine analytically the nu- 
merical obtained exponents. 

The outline of the paper is as follows. In Sec. Ull we 
present both the analytical model to describe sponta- 
neous imbibition and the numerical phase field model. 
In Sec. mil we review the statistical properties of rough 
interfaces describing the anomalous scaling, and in the 
Sec. IIVI we analyze and discuss the obtained numerical 
results. The final conclusions are given in Sec. 



II. THE MODEL 



viscous fluid 

V^v = -Vp, (1) 
A* 

is reduced to the 2-D Darcy's law: 

v{x,y) = -KVp{x,y), (2) 

where v represents the averaged velocity over the di- 
rection z. Since we considerer a constant permeability 
K = 6^/12/i, we can impose the liquid conservation |19l |. 
V • ■»} ~ , obtaining the Laplace equation for the pres- 
sure 



A. The Spontaneous Imbibition phenomenon 

Spontaneous imbibition consists in the advance of a 
viscous fluid of viscosity (like oil, water) in a porous 
medium with air, where motion is driven by capillarity. 
An important classical known result of this kind of mo- 
tion is that the mean advancing front obeys Washburn's 
law [25, H(t) ^ i^/^, thus spontaneous imbibition is a 
case with slowing-down dynamics. 

Experimentally, a porous medium can be well repro- 
duced through the so-called Hele-Shaw cell, which con- 
sists on a two horizontal glass plates of size x L^, 
separated by a small gap b in the 2-direction, where a 
flberglass substrate with copper squares is attached to 
the bottom plate (see Ref. Q) producing a random gap 
spacing. A scheme of the meniscus curvature of the in- 
terface variation due to the quenched noise is depicted in 
Fig. in 
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FIG. 1: Gap variation due to the quenched noise 



The main effect of this random gap is to produce a ran- 
dom capillarity and permeability. However, in our spe- 
cific case of spontaneous imbibition, where the motion 
is controlled solely by capillary forces and the interface 
advances at very low velocity, the fiuctuations due to 
the permeability will be less important than those pro- 
duced due to the random capillarity. Actually, it has 
been shown that the permeability term depends always 
on the interface velocity [la . ll9ll25| . In this sense, we can 
suppose just the random capillarity as the main destabi- 
lizing term. Note that this approximation could not be 
valid in the forced-fiow imbibition case, where the inter- 
face advances at constant velocity and thus, the random 
permeability could become an important term, depend- 
ing on the imposed velocity. 

With this approximation in mind together with the 
fact that h << Lx,Ly, the equation of motion for the 



V^p = 0. (3) 

The boundary conditions 

Ap\i„t = a-K + r]ix,h), (4) 
p{x,y = 0)^Pa, (5) 

complete the equations to describe the spontaneous im- 
bibition phenomenon. Eqns. (O , © , Q and Q are 
the equations of the macroscopic model. Eq. is the 
local thermodynamic equilibrium condition at the inter- 
face, where a is the surface tension of the liquid, k is 
the local curvature of the interface and ri(x, h) represents 
the random capillarity due to fluctuations in gap spacing, 
which it can be approached as rj{x,h) ^ 2acoa9/b. We 
also assume that the contact angle 9 does not vary. 

On the other hand, Eq. ((SJ is the applied pressure 
at the origin of the cell and it is the main difference 
between the spontaneous imbibition and the forced-flow 
imbibition, since in the last case, a constant velocity is 
imposed. 



Washburn's law 

As a simple case of constant capillarity, r/{x,y) = rjo, 
the problem is reduced to compute Eqs. Q and ^ for 
the advance of a flat interface with the boundary condi- 
tions, Eqs. 0) and ©, where now k = 0; obtaining the 
exact solution 

dim ^ KjTJo + Pa) 

dt H{t) ' ^ ' 

where H{t) is the position of the interface. Integrating 
the last expression one obtains the so-called Washburn's 
law 

H{t) = ^Hl + 2at, (7) 

where a ~ Kiijo + Pa) and we have supposed an initial 
height of liquid Hq. Therefore, the expected behavior 
H{t) ~ t^/^ is recovered after an initial transient regime. 
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Linearized Interface Equation 



B. The Phase Field Model 



The general problem of a random field y) is more 
complex to solve. The interface equation can be de- 
rived from Eqs. (O, and 10) using Green func- 
tion analysis. In this case, the Green function G(r, r') 
obeys V^G'(r,r') = 5{r~r') evaluated in the semi-plane 
O = {a;,0 ^ 2/} with Dirichlet boundary conditions. In 
order to get an expression for interface fluctuations, the 
Green equation 

dr' Ur')V'2G(r, r') - G(r, r')V''^p{r') ^ 

L 

[ ds' ■p{s')V'Gis,s') - [ ds' ■Gis,s')V'p{s'), (8) 

Jsl JSl 

is integrated over the volume of the liquid : {a;,0 ^ 
y ^ /i(a;,t)}, where h{x,t) is the interface position, and 
the Darcy equation, Eq. is imposed at the inter- 

face. Linearising in small deviations of the height h{x,t) 
around their mean value H(t) in the Fourier space, one 
can obtain the linearized interface equation 



-aK\q\q% ~ H{t)\q\hg + K\q\fjg, 



(9) 



where we have supposed that correlations grow up slower 
in time than the position of the averaged interface over x, 
which lies in the condition \q\H(t) » 1. The noise term 
in Fourier space is taken as 77, = / dxe~'^''^ri{x, h{x,t)). 
The zero mode of the interface position follows the 
Washburn equation, Eq. ©, in the form H(t) = 
a/H(t), where now, a depends on the noise average 
{i]{x, h{x, t)))x instead of 770. 

The linearized equation, Eg. is the same equa- 

tion obtained in Ref. [20j and contains two mecha- 
nisms to damp the interface fluctuations produced by the 
quenched noise fjq. At small scales, the surface tension a 
controls the fluctuations introducing a dynamic exponent 
z = 3. Actually, it is the regime observed experimentally 
in Ref. 0, indicating that a linear description is valid 
in that situation. The long scales are damped by the 
advancement of the front H{t). These two effects are 
separated by a dynamic crossover length: 
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aK \ 1/2 



H{t) 



(10) 



Since the mean height of the interface follows the Wash- 
burn law, this crossover length increases with tim e, £ x ^ 
t^/^. The numerical results presented in Ref. |2(l . l21| 
shows that this crossover length acts as cut off for the 
fluctuations growth due to the interface is always asymp- 
totically flat on length scales larger than • Thus if the 
velocity of the interface is large enough, the roughening 
of the interface is controlled by which gives rise to a 
dynamic exponent z = 4. 



Our numerical approach is based on the phase field 
model, already introduced in Ref. , where a conserved 
field 4> is used to represent the two phases of the prob- 
lem: liquid/air with the equilibrium values (peq = +1/^1 
respectively. The dynamics of the field is controlled 
by a continuity equation based on a time-dependent 
Ginzburg-Landau model with conserved order parameter 
(model B in 113 ) = VMV^ where = ST/Scf) is 

the chemical potential and the free energy is taken with 
the form / rfr(y(</>) + {eV<j))^/2). A double weU 

potential is also chosen with a linear random term, which 
destabilizes one of the phases and forces the interface to 
advance: 



1 



1 



+ -(/)' - ?/(r)(/). 



(11) 



Here, the noise 77 (r) plays the capillarity role. In our 
numerical integrations we use a dichotomic noise dis- 
tributed in a 2-D system with the values: 



V = 



Vo 



1-»)A 



(12) 



The noise is imposed in the same way that in the experi- 
mental work in a Hele-Shaw cell 0, and it is distributed 
in such way that the sites with 770 occupies the 65% of 
the whole system. Even though in such experiments cap- 
illarity arises as a 3-D effect, wc can relate our noise pa- 
rameters with the Hele-Shaw gap spacing & as r;o ~ 1/6 
and rjA ~ d/b. So the 77^ parameter controls the contrast 
between both noise values whereas the 770 corresponds to 
the capillarity of the cell with a fixed gap. 

The resulting equation for the phase field reads 

= VMV[-(j> + (l)^-e^W^(f>-ij{r)], (13) 

where M is a parameter which we take constant at the 
liquid phase (</> > 0) and zero at the air phase (0 < 0). 
In order to recover the macroscopic results described by 
Eqs. (|2J), Qj O and Q, a matched asymptotic expan- 
sion in the Sharp Interface Limit, e — *■ 0, is necessary. 
This expansion has been done in Ref. [l^. Although in 
that case the phase field is used to describe the Saffman- 
Taylor problem, a similar expansion is also valid in our 
case of a stable interface, with the identification of the 
phase field parameters as: 



where fii is the first order in e term of the chemi- 
cal potential in the asymptotic expansion and — 
—(f)eqtanh{w / \/2) is a kink solution of the phase field 
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equation for a flat interface in equilibrium, and it corre- 
sponds to the zero order in e term of the order parameter. 
The variable w is an inner coordinate of the interface, in- 
troduced in the expansion, and it is perpendicular to it 
at any point. 

Since the chemical potential is equivalent to the pres- 
sure, it has to be fixed at the origin of our system to 
reproduce spontaneous imbibition, = 0) = fXa- In 

Fig. [5]thcrc is plotted a scheme of the chemical potential 
variation along the y-direction. The imbalance created 
in the chemical potential by -t- 77 causes the interface 
to advance following the Washburn law, Eq. iQ. 



1^ 







FIG. 2: Variation of the chemical potential in an arbitrary 
time along the y-direction. The velocity of the interface is 
related to —dyfi. 

Throughout the paper we call applied pressure to the 
value Ha- We will also suppose an initial height of liquid 
Hq to get a low initial velocity, Eq. ©. 



III. SCALING OF ROUGH INTERFACES 

The statistical properties of a one-dimensional inter- 
face defined by a function h(x, t) are usually described in 
terms of the fluctuations of h. The Family- Vicsek scaling 
hypothesis ensures the existence of a lateral correla- 
tion length £c expected to increase as £c t^^^ until it 
reaches the system size ic = L, which it defines a sat- 
uration time tg ■ On the other hand, the vertical 
correlation length or interface global width, defined as 

W{L,t) = {[h{x,t) -h]'^)^/'^, increases as WiL,t) - ff^ 
for t < ts and becomes constant W{L,t) ~ for 
t ^ tg. Here, (..) denotes average over different noise 
realizations and the overbar is an spatial average in the 
X direction. Alternatively, one may study the correla- 
tions over a distance £ <^ L through the local width, 
w{e,t) = {{[h{x,t) - {h)if)e)'^/'^, where (...)£ denotes an 
average over x in windows of size £. Since in the Family- 
Vicsek scenario the local scales grows in the same way as 
the global ones, the scaling behavior of the surface can 
be obtained by looking at the local width. 

Another useful quantity also used to obtain the rough- 
ness exponent a is the power spectrum of the inter- 



face S{k,t) = {hk{t)h-k{t)) , where hk{t) is the interface 
Fourier transform. In the Family- Vicsek assumption it 
scales as 



where s pv is the scaling function 



SFviu) 



const, if u 3> 1 
if u < 1. 



(15) 
(16) 



However, experimental results and several growth models 
have appeared in the last decade showing that global and 
local scales are not equivalent. It is the so-called anoma- 
lous scaling p^ . In this sense, a new scaling ansatz is 
needed to characterize both global and local growth. 

This new scaling ansatz was proposed by Ramasco et 
al |0| and it was a generalization of the Family- Vicscck 
scaling function, Eq. 1)16(1 . 

5(fc,<) = fc-(2"+i)s^(fctiA)^ (17) 
where now the scaling function has the general form 



sa{u) 



if u < 1, 



,2a+l 



(18) 



being ag the spectral roughness exponent. This scaling 
behavior of the power spectrum leads to the following 
dynamic scaling for the local width: 



with the corresponding scaling function 



5(") 



u"'°<= if u < 1 
const, if u ^ 1, 



(19) 



(20) 



where aioc is the local rough exponent and it character- 
izes the roughness at small scales. One of the implications 
of the anomalous scaling is that the local width saturates 
when the system size saturates as well, i.e, at times tg 
and not at the local time ti ~ £^ as occurs in the Family- 
Vicsck scaling. There is an intermediate regime between 
ti and tg where the local width grows as w(£, 
with (3* = j3 - aioc/z. 

Different types of scaling arise from the general scaling 
function, Eq. (|17|l . For < 1 it is always accomplished 
that a/oc = Oig\ in this case the Family- Vicsek scaling 
is achieved when aioc = a, and the intrinsic anomalous 
scaling appears when a > aioc- The main feature of such 
an anomalous scaling is the existence of a temporal shift 
in the power spectrum due to the difference between aioc 
and a. On the other hand, for > 1 always occurs that 
aioc — 1 and the super-rough anomalous scaling appears 
when a = a^. 



IV. NUMERICAL RESULTS 

The numerical study of the imbibition phenomenon is 
made by means of the integration of the phase field equa- 
tion, Eq. (|13|l . in a rectangular lattice of Ax = 1 grid 
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FIG. 3: Squared averaged interface height position, as 
function of time, follows Washburn law for different applied 
pressure fj,a- 



space, with periodic boundary conditions in the x di- 
rection. We have also considered a smaller grid space 
Ax = 0.5 and we reproduce identical results as we ex- 
plain below. The interface position h(x, t) is approached 
by a linear interpolation of the zero of the phase field, 
4>{x, h{x,t);t) = 0. In all of our numerical integrations 
we will suppose an initial height Hq 20 and the noise 
value ?7o = 3.0. Our study is focused to obtain the differ- 
ent scaling exponents and its dependence on the applied 
pressure, /Zq, and on the contrast between noise values 

VA- 



A. Variation in the applied pressure 

Different applied pressures at the origin have been 
studied. In all of them the Washburn law, Eq. |(7)l, is well 
accomplished as we can see from Fig. 01 As a first case 
to analyze, we choose the numerical parameters r]A = 0.3 
and /ia = —3.0. In Fig. 0] there is plotted the temporal 
evolution of the interface. Note that the applied pres- 
sure has the same value in module than the noise value 
770 and thus there are grid sites where the interface is lo- 
cally pinned (Fig. |2J|. In this situation, the interface has 
a slow motion and advances in avalanches. 

Using the scaling concepts described in the last sec- 
tion, we do a statistical description of the interface fluc- 
tuations. Fig. From the interface local width we ob- 
tain the global growth exponent /3 ~ 0.49 ± 0.03 and 
the local one (3* — 0.20 ± 0.04, clearly different from 0, 
indicating the existence of anomalous scaling. On the 
other hand, from the interface power spectrum we ob- 
tain the spectral roughness exponent as = 1.0 ±0.1. The 
temporal shift observed in the power spectrum indicates 
that a > as and since as is not greater than 1, the re- 
lation aioc = tts should be accomplished Therefore, 
the global roughness exponent a can not be achieved di- 
rectly from the power spectrum but it has to be obtained 
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FIG. 4: Temporal evolution of the interface at equal time 
intervals. Numerical parameters are 770 = 3.0, tja = 0.3 and 
Ha = -3.0 



by means of a collapse of the interface local width func- 
tion, Eq. (|19|l . The result of the best collapse is shown 
in Fig. from which we obtain that a = 1.50 ± 0.02 
and z = 3.0 ± 0.1. This plot also shows a local roughness 
exponent of aioc = 0.95 ± 0.03, which is very close to the 
value of the roughness exponent, as- These set of scaling 
exponents are compatible within the intrinsic anomalous 
scaling scenario, and they are clearly different from those 
obtained in Ref. pol |. 

Let us now increase the interface velocity by going up 
the applied pressure, /Xa- Following the same statistical 
procedure as before, we compute the whole set of scaling 
exponents for different values of > —3.0. Fig. [Blshows 
the global width W{L,t) in each case. We observe a 
crossover from (3 ~ 0.5 {fia ~ —3.0) to /3 ~ 0.37 {jia = 
—2.7). In the specific case of /ia = —2.7, Fig. [7| the 
different scaling exponents are a = 1.5 ± 0.02, z = 4.0 ± 
0.1, /? = 0.37 ± 0.03, and aioc = 0.95 ± 0.03. 

Finally, the results obtained in Ref. [23|, a = 1.25, 
z = 4, /? ~ 0.32 and aioc — 1, are recovered when we 
increase the applied pressure up to ~ —1.0. So it is 
important to remark that decreasing the applied pres- 
sure produces three important effects: 1) the roughness 
exponent a changes from 1.25 to 1.5, 2) the interface 
local growth change from a super-rough to an intrinsic 
anomalous scaling description and 3) the dynamic expo- 
nent z = 3 obtained in experimental work is attained 
when the interface advances at low velocities. In Table |l] 
we can see the different scaling exponents computed for 
different values of fia- 

The meaning of large and low velocities has to be un- 
derstood as follows. The interface advances at low ve- 
locity when the pinning effects become important and 
the crossover length is large enough to observe an ini- 
tial regime described by z = 3. Therefore, in this sense, 
we call low velocities the cases of fia < —2.0 and large 
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FIG. 5: Statistical description of interfaces when [i,a = —3.0 
and r\A ~ 3.0: (a) Log-log plot of the interfacial width as 
function of time in different windows sizes £. The data fit 
shows a growth exponent /3 — 0.49 ± 0.03 and the local 
one /3* = 0.20 ± 0.04. (b) Collapse of the interface local 
width function, w{£,t), which is scaled by versus £/t^^°', 
where a = 1.50 ± 0.02, and then z — 3.0 ± 0.1 follows from 
z = a/ p. The picture also shows a local roughness exponent 
aioc = 0.95 ± 0.03. (c) Power spectra of the interface h{x,t) 
at different times. It shows that as — 1.0 ± 0.1. 



FIG. 6; Computed global width W{L,t) by changing the 
applied pressure Ha- The lines show a fit data which allows 
us to obtain the P exponent 
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TABLE I: Scaling exponents obtained when the applied pres- 
sure fia decreases. The contrast of the noise values is kept 
constant rjA ~ 0.3. 



velocities the other ones, > —2.0. 

Concerning the low velocities regime, there is some nu- 
merical inaccuracy in the determination of the interface 
position due to the slow advance of the interface. In or- 
der to check our numerical results, we have performed 
the case of fia — —3.0 with a smaller grid space. We 
have used a Ax = 0.5 and no difference between the case 
of Ax = 1.0 has been observed. 



The crossover and correlation length 

So far we have seen that an initial regime, dominated 
by the surface tension term of Eq. , can be observed 
if the interface velocity is low enough. The point is that 
the correlation length ic = A ■ t^^^ is smaller than the 
crossover length during the time interval studied. In 
order to check it, we are interested in estimating the 
crossover and correlation length. First, we suppose that 
in the case of fia = —2.7, in the regime of z = 4, the 
crossover length, = ^ ■ t^^^, acts effectively as the 
correlation length as Dubc et al have pointed out [2?il |. 
where the constant A can be approached from the inflec- 
tion point of the width collapse. Fig. [7|d, being A ~ 2.5. 
On the other hand we know that the crossover length has 
the expression ^x = C(2/a)^/^ • i^^^, obtained from Eq. 
(|10|l in the limit < 2at, where a follows from Wash- 
burn law, Eq. With these two expression we can 
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FIG. 8: Estimation of the crossover and correlation length 
when the dynamic exponent z = 2> m observed. 
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FIG. 7: Statistical description of interfaces when = —2.7 
and r]A = 3.0: (a) Log-log plot of the interfacial width 
as function of time in different windows sizes i. The data 
fit shows a growth exponent /3 = 0.37 ± 0.03 and the lo- 
cal one /3* = 0.14 ill 0.04. (b) Collapse of the interface local 
width function, w(£,t), which is scaled by versus i/t^^", 
where a = 1.50 ± 0.02, and then z — 4.0 ± 0.1 follows from 
z = a/ p. The picture also shows a local roughness exponent 
oioc = 0.95 ± 0.03. (c) Power spectra of the interface h{x, t) 
at different times. It shows that as — 1.0 ± 0.1. 



obtain the constant C ~ 1.4, which allow us to evaluate 
the crossover length in its full expression, Eq. (|10|) . when 
fia = —3.0. In such case the correlation length grows up 
as £c = 2 • t^/^, where now the parameter 2 has been esti- 
mated from the width collapse in Fig. [SJj. These results 
arc depicted in Fig. |H1 where we can observe that during 
the studied time interval, the condition > -^c is ac- 
complished, and thus the obtained dynamical exponent 
z = 3 is fully justified. 

It should be noted from Eq.© that the crossover 
length takes a certain time t = Hq /2a to reach its asymp- 
totic form — t^^^- It could produce some effect on the 
interface scaling as long as tx < t, where tx is defined 
as the time when ic = £,x- As we can see in Fig. |S1 
in both cases of = —2.9 and fia = —2.7, there is a 
regime, during the transition from z = 3 to z = 4, where 
the interface global width does not increase. Comput- 
ing the parameter a from a fit of the Washburn law, it 
turns out that t ~ 2 • lO'' and 10"^ for fia = —2.9 and 
—2.7 respectively, indicating that t ^ tx in both cases. 
For larger values of the applied pressure this effect occurs 
in an initial stage of the interface growth and it is not 
observed. 

An interesting point to comment from these numer- 
ical results is that the roughness exponent a does not 
change its value when the correlation length catch up the 
crossover length, that is, the growth exponent (3 varies 
accordingly with the change of the dynamic exponent 
from z = 3 to z = 4 keeping constant the a value (see 
Fig. El). Therefore, it lead us to think that the crossover 
length does not change the scaling properties of the in- 
terface except for the dynamic exponent z. It can be 
very useful to extract information from the linear Eq. 
©. Since the initial interface velocity has a finite value 
greater than zero, the condition tx > is always ac- 
complished and therefore there will be always an initial 
regime where fluctuations are dominated by the surface 
tension. Note that this initial regime could be very short 
for a large initial velocity of the interface. Hence, we can 
suppose that the scaling exponents a and a/oc, observed 



when the crossover length is controUing the fluctuations 
growth, are the same than those corresponding to the 
initial regime described by z = 3. 

Intrinsic anomalous scaling in the low velocity regime 

Our numerical results show the influence of the applied 
pressure on the scaling exponents. We have seen that 
the intrinsic anomalous scaling appears when we decrease 
the applied pressure. Actually, the large local slopes, 
which appear when the interface velocity is decreased, 
cause that the scaling properties of the interface follows 
an intrinsic anomalous description. 

In order to approach analytically the global exponents 
in the case of low velocities we should know the scal- 
ing properties of the noise. Likewise, the fact the in- 
terface is locally pinned may introduce some difficulties 
on the quenched noise and we can not know a priori 
how the noise scales. However, the roughness exponent 
a = 1.5 seems to indicate that the interface feels an effec- 
tive white columnar noise, that is ?7(x) with correlations 
{r]{x)rj{x')) ~ 5{x — x'). Indeed, if we rescale the lin- 
earized interface equation, Eq. by 

x^hx h^b'^h bH, (21) 

the noise scales then as 

Tj{x) ^ r^77(x), (22) 

and the roughness exponent we obtain is a ~ 1.5. This 
roughness exponent is the same value obtained in the 
local QEW model with columnar noise [s^l • 

It should be pointed out that this assumption of colum- 
nar noise is related to the number of sites where the in- 
terface is pinned; note that the fractions of the interface 
which are locally pinned increases as the applied pressure 
decreases up to values lower than fia < —3.0 where the 
interface does not advance at all. So in this sense, the 
roughness exponent a = 1.5 is observed when the inter- 
face advances near a critical situation of pinning. These 
results agree with the hypothesis of Parisi (2^ who sug- 
gest that near the pinning transition, both QEW and 
QKPZ models follows the same universality class than 
in the QEW model with a columnar noise r]{x). How- 
ever, we have to remark that only the roughness expo- 
nent a = 1.5 obtained in our model adjusts with such an 
universality class. 

On the other hand, when the applied pressure is in- 
creased and therefore we move further away from the 
critical situation we recover the results of a = 1.25 corre- 
sponding to the university class of the local QEW model 
with quenched noise ?7(a;,^)0]. In this sense, the other 
observed roughness exponent a ~ 1.34 can be understood 
as an intermediate regime between the value a = 1.25 
and a = 1.5 when the applied pressure is decreased. 
Likewise, in the next section we will discuss an analytical 
approach which describes the quenched noise scaling in 




FIG. 9: Statistical description of interfaces described by Eq. 
11311 with a gaussian noise distributed at the 100% of the sys- 
tem with mean (rj) = 3.0 and strength A77 = 0.6. The applied 
pressure is fia ~ —2. 5. (a) Collapse of the interface local width 
function, w{l,t), which is scaled by t" versus ^/t"/", where 
a = 1.51 ± 0.02, (3 = 0.51 ± 0.03 and then z = 3.0 ± 0.1 fol- 
lows from z = a//3. The picture also shows a local roughness 
exponent aioc = 0.9 ± 0.03. (c) Power spectra of the interface 
h{x,t) at different times. It shows that Qs = 1.0 ± 0.1. 



the high velocity regime, when the values a ~ 1.25 and 
a = 1.33 are observed. 

In order to check the validity of this numerical results, 
we have performed the simulations with another kind of 
noise. We have used a gaussian field distributed at the 
100% of the system with mean (77) = 3.0 and strength 
A77 = 0.6. When the apphed pressure is /Iq = —2.5, the 
interface advance at low velocities and we recover the 
same results of a = 1.5 and z ~ i within the intrinsic 
anomalous scaling, Fig|^ Therefore, the columnar noise 
assumption is also correct when we use a different noise 
distribution. 

Comparing to the experimental results, our a = 1.5 
differs from the value a = 2.0 obtained in Ref. 0. How- 
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FIG. 10: Effect of increasing the contrast between the two values of the noise in the interface profile. The values used are: (a) 
rjA = 0.3 and (b) rjA ~ 0.5. The main difference is the switch from super-rough scaling to intrinsic anomalous one. 



ever, the other independent exponents, aioc and z are 
in good agreement. This discrepancy in the a exponent 
could be attributed to a very low velocity of the exper- 
imental interface, not attained in the numerical integra- 
tion of the phase field, introducing a different propcrtie 
on the scaling of the quenched noise; however, it remains 
as an open question. 



B. Variation in the capillary contrast 

Another interesting issue to study is what happens 
when we increase the contrast between the two values 
of the noise keeping constant the applied pressure. Actu- 
ally, it is similar to increase the strongness of the capil- 
larity of the Hele-Shaw cell and a similar study was car- 
ried out experimentally by Soriano et al 0. So in this 
sense, we are interested to study the change of the scal- 
ing exponents when we increase the parameter rjA- Our 
starting point is the known results of super-rough inter- 
faces with a = 1.25 and z = 4 obtained in the case of 
Ho = —1.0 and ija = 0.3. In Fig. ^|we can see two differ- 
ent interface profile with two different capillary contrast, 
•qA = 0.3 and tja = 0.5 respectively, where we observe the 
morphological difference between both interfaces. While 
in the former case the interfaces are characterized by a 
smooth profile, large local slopes appear when the cap- 
illary contrast is increased. Note that due to numerical 
limitations, the capillary contrast is not able to be too 
large. Actually, for values -qA > 0.6 we have observed 
that the interface becomes unstable and a clear statis- 
tical description of the interface fluctuations can not be 
performed. 



The different scaling exponents for the case of rjA = 0.5 
are computed in Fig. 1111 obtaining: a = 1.33 ± 0.03, 
z = 4.0±0.1, P = 0.33±0.03, (3* = 0.16±0.04 and aioc = 
0.70 ± 0.03, which are compatible within the intrinsic 
anomalous scaling framework. In Table ^] we can realize 
how the scaling exponents change by varying the noise 
contrast. While the global exponents (a, (3 and z) are 
quite similar in both cases, the interface local behavior 
changes from a super-rough to an intrinsic anomalous 
description. 

If we use a gaussian noise instead of a dichotomic noise 
as in the previous section, we obtain a similar variation in 
the local scaling properties of the interface fluctuations. 
It turns out that increasing the strength of the gaussian 
noise, Arj, an intrinsic anomalous scaling of the interface 
fluctuations is observed. Therefore, there seems to be a 
relation between the capillary contrast of the dichotomic 
noise and the strength of the gaussian field. Likewise, 
Laurila et al |23| have observed a slightly dependence of 
the strength of gaussian noise on the roughness exponent 
a in a similar way that our results of Table ^ when we 
change the capillary contrast of the dichotomic noise. 
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TABLE II: Scaling exponents for two different capillary con- 
trast, rjA, with an applied pressure fia = —1.0. 
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FIG. 11: Statistical description of interface with an increased 
capillary contrast, r]A ~ 0.5: (a) Log-log plot of the interfacial 
width as function of time in different windows sizes £. The 
data fit shows a growth exponent /3 = 0.33 ±0.03 and and the 
local one /3* = 0.16 ± 0.04.(b) Collapse of the interface local 
width function, ■w{£,t), which is scaled by versus £/t^^°', 
where a = 1.33 ± 0.03, and then z = 4.0 ± 0.1 follows from 
z — a/ p. The picture also shows a local roughness exponent 
otioc = 0.70 ± 0.03. (c) Power spectra of the interface h{x, t) 
at different times. It shows that as — 0.75 ± 0.05. 
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FIG. 12: Dichotomic noise, r]{x,h), evaluated in a fixed in- 
terface h{x,tf) for each x in the case of the capillary contrast 
riA = 0.5. Ac represents an averaged length where the noise 
takes a constant value. 



Scaling properties of the quenched noise 

Our objective now is to try to determine analytically 
the global scaling exponents observed in each case of the 
studied capillary contrast. While the dynamic exponent 
z can be well understood from the different terms of the 
linear equation, Eq. (PJ, the roughness exponent a in- 
volves the knowledge of the scaling of the noise and it 
can be a more difficult task. 

First, we assume the general form for the noise corre- 
lation 



{rjix,h)rjix',h')) c^f{x,x')gih,h') 



(23) 



In order to describe the scaling of the noise, we are inter- 
ested in find out the expression for the functions /(x, x') 
and g{x,x'). To do it, we first focus on the case of 
7]A = 0.5 (Fig. 10(b) I. In this situation, an important 
feature of the interface velocity is that it has a great in- 
crease when the interface cross through the sites where 
the noise takes the largest value of Eq. (|12|l . In some 
way, interface seems to avoid these sites and it prefers to 
be where the noise is lower, that is r\ = 770. It leads to 
think that there has to be a length Ac below which the 
noise has a constant value. Since there is no imposed ve- 
locity but the interface moves only due to the capillarity 
of the media, this length Ac(i, r\A) is expected to increase 
with time and with the capillary contrast. An interpre- 
tation of Ac can be seen in Fig. ^1 where the values of 
the quenched noise r\{x, h) in a fixed interface h(x, tf) are 
plotted, being an arbitrary chosen time. Note that Ac 
has to be understood as an average length in such way 
that we can suppose the function f(x,x') ~ const for 
\x — x'\ < Ac and f{x,x') ~ 5{x — x') for \x — x'\ > Ac. 

Since the only length scale characterizing the fluctua- 
tions growth in the x direction is the correlation length 
£c, two different scaling for the noise emerge depending if 
£c > Ac or £c < Xc- In Fig. E| there is plotted the length 
Ac in the two cases of capillary contrast, rjA = 0.3 and 
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0.5 is the same that the value obtained in the Sec. IIV Al 
in the case of fXa = —2.0 (see TablePl, indicating that, in 
that case, the scaling of the noise has to follows also Eq. 
(|24|l . Therefore, we analytically obtain the roughness 
exponents in the high velocity regime when we vary either 
the capillary contrast or the applied pressure. 

Eventually, since the crossover length controls the in- 
terface fluctuations in both cases, we conclude that the 
dynamic exponent has to be z = 4, and then we achieve 
the whole set of global scaling exponents described in 
Table ini 



FIG. 13: The length Ac of the quenched noise is always greater 
than interface correlation length Ic in both cases of 77A = 0.3 
and 0.5. 

rjA = 0.5, compared with the interface correlation length 
£c — t^^"^. We can see that in both cases Ac > (-c and 
thus we can take the correlation of the noise along the 
x-direction as a constant. 

Concerning the expression for the function g{h^ h'), we 
will suppose the simplest case of g{h,h') — 5{h — h'). 
However, note that two different kind of scaling emerge 
depending if the quenched noise is evaluated on the in- 
terface T][x, h) or on the mean value of the interface, 
r]{x,h) ~ r]{x,H{t)). Therefore, we conclude that the 
quenched noise has to scale as 

r]{x,h) ^b°b-'iT]{x,h), (24) 
ri{x,H) ^b°b~^r]{x,H), (25) 

when we rescale it by Eq. (|21|l . The term with 5° comes 
from the constant function f{x,x'). We have also as- 
sumed the Washburn law H{t) ~ t^/^. 

It should be noted that these noise scalings are not 
valid when the interface movement is near a critical pin- 
ning situation ( as in the previous section where /ia ^ 
—3.0) since in such a case there is an important number 
of sites where the interface is locally pinned and the ef- 
fects of a columnar noise, Eq. H22|) . become important. 



Global scaling exponents 

Once we know how the noise scales, we can compute 
analytically the global scaling exponents in the studied 
cases of capillary contrast, tja = 0.3 and ija = 0.5. We do 
a scaling analysis of the linear equation, Eq. taking 
into account only the surface tension term. Imposing 
the noise scaling as Eq. 1)24(1 . we obtain the exponent 
a = 1.33, which is in accordance to the case of tja = 0.5. 
On the other hand, if we suppose that noise scales as Eq. 
(|25|l . the exponent is then a = 1.25, which adjusts the 
numerical value obtained in Ref. |20j | and in our case of 
T]A = 0.3. Note that the a = 1.33 observed when tja = 



V. CONCLUSIONS 

We have presented numerical results of spontaneous 
imbibition using a phase field model. Our study has been 
focused on the dependence of the scaling exponents when 
we vary either the applied pressure imposed at the origin 
of the system or the capillary contrast of the dichotomic 
quenched noise. 

By decreasing the applied pressure, we have got in- 
terfaces advancing at low velocities. In this situation of 
spontaneous imbibition, our results are in good agree- 
ment with the experimental work made in Ref. 0. In- 
deed, our numerical study confirms that there is a regime 
of low velocities where the growth of the interface fluctu- 
ations is controlled by the surface tension of the interface, 
following £c — t^^^. We have also observed that the cut 
off crossover length, , acts as an effective correlation 
length when the initial interface velocity is higher, ac- 
cordingly to the results of Dubc et al |23 . 

On the other hand, our results show that the local 
scaling properties of the interface change as we decrease 
the applied pressure. More precisely, we obtain an in- 
trinsic anomalous description at the regime of low veloc- 
ities, where the interface is locally pinned. In addition, 
the roughness exponent we observe in such a situation 
is a ~ 1.5 which is the same value obtained in the local 
QEW model with a columnar noise, ri{x). 

It is also important to remark that the intrinsic anoma- 
lous scaling is not related to the low velocity of the in- 
terface but to the capillary contrast of the noise. Indeed, 
when we increase the contrast between both noise values 
keeping constant the applied pressure, the local fluctua- 
tions description of the interface changes from a super- 
rough to an intrinsic anomalous scaling, indicating the 
existence of a direct relation of the capillary contrast and 
the intrinsic anomalous scaling, which had been observed 
experimentally in Ref. 0| . We have also introduced some 
scaling properties of the dichotomic quenched noise which 
have allowed us to determine analytically the obtained 
numerical scaling exponents. 
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